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Tltis  letter  reports  on  a  method  of  maximizing  the  power  output  from  highly  saturated  cladding-pumped  continuous-wave 
single-frequency  fiber  amplifiers  simultaneously  taking  into  account  the  stimulated  Brillouin  scattering  and  transverse  modal 
instability  thresholds.  This  results  in  a  design  figure  of  merit  depending  on  the  fundamental  mode  overlap  with  the  doping 
profile,  the  pea!:  Brillouin  gain  coefficient,  and  the  peak  mode  coupling  gain  coefficient.  This  figure  of  merit  is  then 
numerically  analyzed  for  three  candidate  fiber  designs  including  standard,  segmented  acoustically  tailored,  and  micro 
segmented  acoustically  tailored  photonic  crystal  libers.  It.  is  found  that  cadi  of  the  latter  two  fibers  should  enable  a  50% 
higher  output  power  titan  standard  photonic  crystal  fiber,  ©  201 4  Optical  Society  of  America 
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High  power  continuous-wave  single  frequency  sources  are 
important  for  metrological  applications  such  as 
gravitational  interferometry  [l],  coherent  laser  radar  [2], 
and  non’ linear  inversion  [3].  The  maximum  power 
achievable  with  single  frequency  fiber  amplifiers  is 
limited  either  by  stimulated  Brillouin  scattering  (SBS)  or 
thermal  modal  instability.  Recently,  a  single  frequency 
fiber  amplifier  with  an  SRS-fimitcd  output  power  of  811 
Watts  was  demonstrated  using  a  photonic  crystal  fiber 
with  gain  and  segmented  acoustic  tailoring  (SAT)  [4].  The 
gain  tailoring  of  this  fiber  primarily  served  to  mitigate 
transverse  modal  instability  while  the  acoustic  tailoring, 
together  with  a  longitudinal  thermal  gradient,  served  to 
mitigate  SBS,  The  acoustic  tailoring  was  devised  in  order 
to  distribute  the  mode  field  intensity  approximately 
evenly  over  regions  of  differing  acoustic  velocities  in  order 
to  reduce  the  effective  peak  Brillouin  gain  coefficient.  It 
has  been  shown  theoretically  that  there  is  an  optimal 
acoustic  domain  size  for  SBS  suppression  in  thermally 
uniform  fibers  that  exploits  elastic  coupling  between 
reduced-size  acoustic  domains  to  suppress  the  resonant 
electrostrictive  response  insulting  in  a  smaller  peak 
Brillouin  gain  coefficient  [5]. 

A  sizeable  amount  of  evidence  has  been  presented  that 
suggests  the  fundamental  mechanism  leading  to  modal 
instability  is  Stimulated  Thermal  Rayleigh  Scattering 
(3TRS)  [6].  The  SBS  and  8TRS  processes  share  several 
key  aspects  in  common.  They  both  can  be  described  as 
the  scattering  of  one  mode  into  another  due  to  a  non¬ 
linear  refractive  index  response  to  the  intensity 
interference  pattern  of  the  tw  o  modes.  In  the  case  of  SBS 
the  inodes  are  separated  from  each  other  the  -16  GHz 
Stokes  shift  and  power  may  be  scattered  between  the 
same  or  different  transverse  modes  while  in  the  case  of 
STRS  the  modes  are  separated  by  a  frequency  offset 
typically  of  the  order  of  kHz  and  the  power  is  scattered 
between  different  transverse  modes.  In  the  case  of  SBS 


the  index  disturbance  is  caused  by  electrostriction  while 
for  STRS  it  is  caused  by  the  thermo- optic  effect.  For  most 
non-linear  processes,  increasing  the  mode  field  area 
increases  the  threshold.  The  STRS  the  threshold  on  the 
other  hand  decreases  with  increasing  mode  field  area  in 
most  cases  as  will  be  shown  below.  This  sets  up  a 
situation  wherein  design  principles  for  Increasing  the 
STRS  threshold  oppose  those  for  increasing  the  threshold 
of  SBS  as  well  as  most  other  parasitic  nonlinearities.  The 
maximum  combined  threshold  for  single -frequency 
systems  such  as  those  considered  here  may  then  be 
achieved  at  the  design  point  where  the  SBS  and  STRS 
thresholds  coincide. 

This  letter  begins  with  a  derivation  of  the  conditions  for 
maximum  output  power  from  a  single J  freque  ncy  fiber 
amplifier  subject  to  SBS  and  STRS,  followed  by  analysis 
of  example  fiber  designs,  and  the  conclusion. 

The  SBS  threshold  power  is  given  by  [7] 

to 

cff.SBS 

where  Ac(C  is  the  non-linear  effective  mode  area,  gs  is  the 
Brillouin  gain  coefficient,  ZcOBS  is  the  effective  fiber  length 
for  SBS,  and  CSBS  is  a  dimensionless  constant  of 
approximate  value  10-20  depending  on  seeding  conditions 
for  the  backward-propagating  Stokes  wave  and  the 
definition  of  the  SBS  threshold. 

In  order  for  (l)  to  be  useful  it  is  necessary  to  impose  a 
condition  on  total  pump  absorption  to  ensure  adequate 
optical  efficiency 

C^^riN,L^SBS~  ® 

where  A,  s  Ni<mcrap  for  a  highly  saturated  amplifier  and 
where  Mc|ad  is  the  pump  cladding  area,  is  the  doped 


core  area.  f1  is  the  overlap  fraction  of  the  fundamental 
mode  with  the  doped  core  area, 

ri  sJJk01v'2'T  (3) 


and  Nim  is  the  doping  concentration  which  leads  to  the 
following  value  for  the  SBS  threshold 
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Thanks  to  a  recently  published  expression  for  the  modal 
power  transferred  in  an  amplifier  due  to  STRS  [$],  an 
expression  may  be  derived  in  a  similar  fashion  for  the 
STRS  threshold  in  a  well -saturated  amplifier.  The  power 
transferred  through  the  STRS  process  into  a  higher  order 
mode  at  a  single  shifted  frequency  is  given  by 


C(0  =  Co(Oexp 


x4^f\p'(z)dz 


(5) 


where  P2  0  ( L)  is  the  output  power  in  the  higher  order 
mode  in  the  absence  of  inter  modal  coupling,  Xo the 
STRS  mode  coupling  coefficient.  [7],  g  is  the  effective 
saturated  signal  gain,  is  the  small  signal  gain,  and 
P,  (r)  is  the  power  in  the  fundamental  mode.  Borrowing 
the  methodology  from  the  SBS  case,  the  threshold 
condition  is  reached  when  the  term  in  the  exponential 
reaches  a  suitably  high  value.  Counter-pumped 
amplifiers  exhibit  more  uniform  saturation  profiles  than 
ar pumped  amplifiers  due  to  the  similar  behavior  with 
length  of  the  pump  and  signal  intensities.  In  the  case  of  a 
highly- saturated  counter-pumped  amplifier  the 
approximation 


is  very  useful.  This  approximation  is  fairly  accurate  if  the 
residua]  pump  power  is  comparable  to  the  seed  power,  a 
condition  not  uncommon  experimentally.  The  saturated 
gain  is  given  by 


g  =  N- 

O-  101 


or  cr 
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from  which  the  highly -saturated  and  small- signal  limits 
may  be  derived  leading  to  the  following  approximate 
expression  for  the  STRS  threshold  power 
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where  is  the  pump  cladding  area,  ,4d  is  the  doped  corn 
area,  F,  is  the  overlap  fraction  of  the  fundamental  mode 
with  the  doped  core  area,  Nl0ii  is  the  doping  concentration, 
Xo  is  the  STRS  inode  coupling  coefficient,  ZcfmfL,  is  the 
STRS  effective  length  (which  may  differ  slightly  from  the 
SBS  effective  length),  Csm  is  the  STRS  dimensionless 
constant,  and  d  is  a  combination  of  cross-sections  given 
by 


<j  = 


+v.X 


(91 


where  the  subscripts  e  and  a  denote  emission  and 
absorption  while  j  and  p  denote  signed  and  pump.  The 
threshold  coefficient  Qm  can  be  expected  to  differ 
depending  on  the  particular  characteristics  of  the 
amplifier  under  consideration  nevertheless  the  threshold 
dependencies  given  in  (8)  should  be  useful  for  optimizing 
fiber  amplifier  designs  and  configurations.  Once  again 
imposing  the  pump  absorption  constraint  and  using  the 
approximation  Ad  =  F,^.  leads  to 
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(10} 


Comparing  (4)  and  (10)  it  is  immediately  evident  that 
increasing  the  mode  field  area  increases  the  SBS 
threshold  but  decreases  the  STRS  threshold.  If  modal 
area  is  chosen  as  the  independent  fiber  design  parameter, 
then  setting  the  STRS  and  SBS  thresholds  equal  yields  an 
optimal  cladding  area  and  maximum  attainable  output 
power.  The  optimal  cladding  area  is  given  by 
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It  is  very  interesting  to  note  that  the  optimal  cladding 
area  depends  on  the  SBS  or  STRS‘Suppressing  properties 
of  the  core  design  as  reflected  in  the  values  of  gH  and  j0  . 
Substitution  of  (ll)  into  either  (4)  or  ( 10)  yields  the 
maximum  power. 
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(12) 


Thus  the  prescription  for  maximum  power  is  to  make  the 
mode  field  area  as  large  as  possible,  and  then  set  the 
cladding  size  to  the  optimum  value.  The  dependencies  in 
(12)  are  very  intuitive.  For  given  material  parameters, 
the  maximum  power  may  be  increased  by  increasing  the 
fundamental  mode  overlap  fraction,  decreasing  the 
maximum  value  of  the  Briilouin  gain  coefficient,  and 
decreasing  the  maximum  value  of  the  mode  coupling 
coefficient.  Equation  (12)  is  not  expected  to  yield  precise 
values  due  to  the  uncertainties  present  in  the  numerous 
parameters  therein.  Furthermore,  the  experimental 
determination  of  these  parameters  is  beyond  the  scope  of 
this  work.  Rather,  this  expression  may  be  used  to 
compare  candidate  fiber  designs  with  tailored  doping  and 
acoustic  profiles  with  respect  to  power  scaling  potential. 

One  successfully  demonstrated  method  of  profile 
tailoring  is  the  microstructuring  of  the  core  whereby  the 
core  design  is  constructed  of  pixels  of  different  material 
compositions  in  a  prescribed  geometric  arrangement  [4] . 
A  given  set  of  finite  materials  constrains  the  power 
optimization  using  (12).  For  the  purposes  of  comparing 
waveguide  design  approaches,  it  is  convenient  to  use  the 
material  set  previously  employed  [4]  consisting  of  pure 


fused  silica  and  two  varieties  of  aluminum,  germanium, 
fluorine  and  ytterbium- doped  fused  silica.  The  present 
aim  is  then  to  assess  the  power  scaling  potential  of 
different  arrangements  of  pixels  in  the  core  using  this 
material  system  while  keeping  the  other  design 
parameters  fixed.  Under  these  constraints,  the  proper 
figure  of  merit  derived  from  (12)  is 

^  =  U  / -Js/iZ,  ■  (is) 

Vectorial  finite  element  techniques  previously  described 
have  been  used  to  calculate  the  fundamental  mode 
overlap  fraction  and  the  peak  Brillouin  gain  coefficient  for 
the  cases  treated  here  [5,9].  Additionally,  these  have  now 
been  adapted  to  calculate  the  modal  coupling  constant  %0 
as  will  now  he  described. 

Assuming  harmonic  time  dependence,  the  frequency 
domain  heat  equation  governing  the  temperature  within 
the  fiber  is 

iQpCT  (r± ,  n)  =  tf71T(r±>  fl)  +Q{r±,Q)  (14} 

as  has  been  previously  described  [6,8]  where  p  is  the 
mass  density,  C  is  the  heat  capacity,  k  is  the  thermal 
conductivity,  T  is  the  temperature,  0  is  the  heating  term, 
and  Q  is  frequency  of  a  given  harmonic  component  of  the 
oscillating  temperature.  Invoking  the  finite  element 
technique  (14)  begins  the  variational  form  governing  the 
temperature  behavior 


3o  =  jJ  ^%(A,n)[2 +^|vr(r1,£i) 


-e(rltn)r(?i(n) 


(15) 


d2r. 


The  temperature  and  heat  distributions  are  then 
approximated  by  their  values  at  the  vertices  of  the  same 
curvilinear  triangular  mesh  used  in  the  optical 
calculations  and  the  area  integral  approximated  as  a  sum 
of  integrals  over  elements  with  locally  uniform  material 
properties  leading  to  a  linear  system  determining  the 
frequency  dependent  temperature  distribution 

{ iQpCM  +  rit)  Ta  =  LnATfi  =  Q,, .  (16) 


Upon  solving  (16)  for  the  temperature  vector,  the 
temperatures  and  mode  field  profiles  may  then  be 
numerically  integrated  over  the  fiber  cross- section  to  yield 
the  frequency  “dependent  modal  coupling  coefficient 


Zo  (°)  = 


dT  Kn 


i=l 


(17) 


where  dn/dT  is  the  thermo-optic  coefficient,  (1  -;/)  is  the 
quantum  defect  fraction,  ktt  is  the  vacuum  wave  vector, 
and  is  the  fundamental  modal  index. 

Using  (13),  three  candidate  fiber  designs  with  different 
doping  and  acoustic  profiles  are  assessed.  The  first  of 
these  has  uniform  doping  throughout  most  of  the  core. 
The  second  is  similar  to  the  SAT  fiber  described 
previously  [4],  and  the  third  uses  a  segmented 
acoustically  tailored  core  but  with  smaller  acoustic 
domains  that  is  termed  a  micro"SAT  design.  In  order  to 


subdivide  the  core  cells  into  smaller  domains  in  the  micro’ 
SAT  fiber  design  it  is  necessary  to  introduce  a  sLight 
asymmetry  into  the  lattice.  For  consistency  of 
comparison,  this  arrangement  is  also  maintained  ibr  the 
PCF  and  SAT  fiber  designs.  All  optical  and  mechanical 
properties  of  the  fibers  ore  the  same  as  those  given  in 
prior  works  [5,10]  except  that  the  two  varieties  of  Ylr 
doped  glass  are  characterized  by  acoustic  velocities  of 
5851  m/s  and  5769  m/s  and  bulk  Brillouin  line  widths  of 
45  MHz  and  46  MHz  but  have  the  same  Yb  doping 
concentrations. 
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Fig,  1.  Doping  (top  row)  and  acoustic  profiles 
(bottom  row)  of  a  PCF  (left),  SAT  (center)  and  pSAT 
(right)  fibers.  Darker  shades  of  grey  represent 
lower  acoustic  velocities  whereas  all  doping 
concentrations  are  constant. 

Fig.  1  depicts  the  acoustic  and  doping  profiles  of  each 
fiber.  These  include  the  first  baseline  fiber  design 
considered  is  a  photonic  crystal  fiber  (PCF)  with  an 
arrangement  of  seven  uniform  cells  of  doped  material  in 
the  core  taking  the  place  of  air  holes  in  the  photonic 
crystal  lattice  resulting  in  a  doped  area  of  604  square 
micrometers.  The  lattice  pitch  is  1,0  micrometers  with  a 
hole  diameter  to  pitch  ratio  of  0.17  resulting  in  a  mode 
field  diameter  of  27  micrometers  and  a  mode  field  area  of 
576  square  micrometers. 

The  second  fiber  design  considered  is  a  SAT  PCF  with  a 
design  similar  to  one  previously  discussed  [4] .  This  fiber 
has  a  center  cel!  surrounded  by  alternating  cells  of  lower 
velocity  doped  fused  silica  and  urr doped  fused  silica. 
Each  cell  has  an  area  of  86  square  micrometers. 

The  third  fiber  design  considered  is  a  pSAT  PCF  with  a 
new  design.  The  acoustic  domain  size  in  this  case  is  2. 1 
micrometers,  which  is  slightly  larger  than  the  optimal  size 
determined  previously  for  a  randomly  acoustically 
microstmctured  PCF  [&].  In  contrast  to  the  random 
acoustic  structure,  the  uSAT  PCF  has  the  domains 
arranged  in  a  specified  fashion.  Each  core  cell  is 
subdivided  into  19  subcells.  Within  each  cell  there  are  6 
high* velocity  subcells,  4  low- velocity  subcefls,  and  9  un~ 
doped  subcells.  This  arrangement  was  determined 
through  trial  and  error  to  yield  the  lowest  predicted  peak 
Brillouin  gain  coefficient.  It  is  important  to  note  that  for 
acoustic  domains  this  small,  vectorial  analysis  is  required 
to  predict  the  Brillouin  gain  spectrum  to]. 

Figure  2  shows  the  calculated  Brillouin  gain  spectra  for 
the  three  fibers.  The  PCF  exhibits  a  single  large  peak 


near  !  5,94  GHz  with  a  small  residual  peak  at  16*28  GHz 
due  to  residual  modal  overlap  with  undoped  regions  at  the 
core  edge.  The  SAT  PCF  exhibits  three  we 0" defined 
peaks  corresponding  to  the  three  acoustic  velocities 
comprising  the  core  at  15.73  GHz,  15.95  GHz,  and  16.28 
GHz.  These  three  main  peaks  are  also  evident  in  the 
spectrum  of  the  LuSAT  PCF,  but  they  are  slightly 
broadened  and  the  two  lower  frequency  peaks  are 
up  shifted  in  frequency* 

Figure  3  shows  the  calculated  Stimulated  Thermal 
Rayleigh  Scattering  mode  coupling  gain  spectra  for  the 
three  fibers.  The  peak  mode  coupling  gain  occurs  at  a 
frequency  shift  of  -3.1  kHz  that  is  the  same  for  all  three 
fibers.  The  mode  coupling  gain  profiles  are  nearly 
identical  for  the  SAT  and  a  SAT  PCF. 


Fig.  2.  Calculated  Brillouin  gain  spectra  of  the  three  libers 

Table  1  presents  the  relevant  amplifier  performance 
parameters  under  the  conditions  described  above.  Both 
the  SAT  PCF  and  jiSAT  PGF  exhibit  an  improvement  of 
about  50%  in  maximum  power  over  the  PCF* 
Interestingly  even  though  the  pSAT  PCF  has  an  identical 
modal  coupling  constant  to  the  SAT  PCF  as  well  as  a 
lower  peak  Brillouin  gain  coefficient,  its  lower  overlap 
fraction  lowers  its  figure  of  merit  below  that  of  the  SAT 
PCF*  This  shows  that  it  is  important  to  consider  the 
overlap  fraction  of  the  fundamental  mode  with  the  doped 
core  area  when  assessing  the  power  scaling  potential  of  a 
transverse  waveguide  design. 


Fig.  3.  Calculated  mode  coupling  gain  for  the  three  fibers 
where  f  -  2nD. 


Table  1*  Amplifier  Feribrnanoe  Parameters 


Fiber  Type 

F] 

gtilm/W] 

^oll/m] 

F  [W12] 

PCF 

0.90 

1.5*10 11 

0.019 

1,7x10s 

SAT  PCP 

0.57 

4.9x10  i- 

0*010 

2.6xiOG 

fiSAT  PCP 

0*47 

3.5x10- 12 

0.910 

2.5x10s 

The  results  presented  here  show  that  maximizing  the 
output  power  of  single 'frequency  continuous  wave  fiber 
amplifiers  requires  considering  simultaneously  both  the 
thermally-induced  mode  coupling  gain  as  well  as  the 
Brillouin  gain  coefficient.  Furthermore,  the  line  of 
development  presented  here  is  applicable  to  other  non- 
linearities  that  have  thresholds  that  scale  with  the  mode 
field  area  as  the  SBS  threshold  does*  Due  to  gain 
saturation,  the  SIRS  threshold  exhibits  distinctly 
different  behavior,  notably  a  decrease  as  mode  field  area 
is  increased  for  a  given  pump  cladding  size. 
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